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SUMMARY 


The closed-loop stability of linear time-invariant multivariable systems con- 
trolled by linear quadratic (LQ) full state feedback regulators is investigated when 
certain types of nonlinear gains are present in the feedback loop. The nonlineari- 
ties N(a) considered are assumed to violate the LQ regulator stability condition, 
a N(cr) > 0.5 a 2 / either (1) for values of a away from the origin (a = 0) or (2) for 
values of a in a bounded region containing the origin. It is proved that there 
exists a region of attraction for case (1) and a region of ultimate boundedness for 
case (2), and expressions are derived for these regions. The analytical results 
obtained provide methods of selecting the performance function parameters in order to 
design LQ regulators with better tolerance to nonlinearities. The results are 
demonstrated by application to control of the pitch-axis attitude and elastic motion 
of a large, flexible space antenna in the presence of saturation and hysteresis non- 
linearities in the actuators. 


INTRODUCTION 

Synthesis of control systems for state or output regulation of linear multivari- 
able systems is one of the most important problems of system theory. In particular, 
for linear time-invariant (LTI) systems, the infinite duration linear quadratic (LQ) 
regulator consists of constant linear feedback of the system state vector. Under 
certain stabilizability and detectability conditions, the closed-loop system is 
asymptotically stable (ref. 1). The infinite duration optimal LQ regulator has been 
shown in references 1 and 2 to also have highly desirable robustness properties, 
namely, an infinite gain margin, a phase margin of ±60°, and tolerance to single- 
valued, memory less nonlinearities in the feedback loop which belong to the (0.5,®) 
sector. (A function N(a) is said to belong to the (K^,^) sector if N(0) = 0 and 
K < a N(a) < K2 <j 2 for a t 0.) In practice, however, many nonlinearities do not 
satisfy this sector condition. Figure 1 shows some commonly encountered nonlineari- 
ties which violate the (0.5,®) sector condition. For example, a saturating amplifier 
has effective gain of less than 0.5 for a < and a > 1 2 (i« e »/ "away from the 

origin"). For the dead zone shown in figure 1(b), the sector condition is violated 
in a neighborhood of the origin. A hysteresis nonlinearity (with a gain greater than 
0.5 in the linear region) also violates the sector condition in the neighborhood of 
the origin (fig. 1(c)). (Also, since hysteresis has memory, the LQ robustness 
property mentioned previously does not apply.) Therefore, this paper extends LQ 
regulator robustness to a broader class of nonlinearities, namely, (1) to those which 
satisfy the sector condition at least in a neighborhood of the origin and (2) to 
those which satisfy it in regions away from the origin, but perhaps violate it in a 
neighborhood of the origin. 

The closed-loop stability of linear systems with nonlinearities in the feedback 
loop (that is, the so-called Lur’e problem) has received considerable attention in 
the literature. (See ref. 3 for a brief history of the problem.) The nonlinearities 
considered in much of the literature were confined to a particular sector, called the 
stability sector. The stability of systems with nonlinearities that escape the sta- 
bility sector were investigated in references 4 and 5 for single-input, single-output 
(SISO) systems. For SISO Lur ' e-type systems, it was proved in reference 4 that when 



the nonlinearity violates the stability sector (i.e., Popov sector (ref. 3)) away 
from the origin, there exists a region of attraction. That is, all the trajectories 
originating in that region asymptotically approach the origin. A method for obtain- 
ing an estimate of the regions of attraction was also given in reference 4. The 
analysis was extended to multi-input, multi-output (MIMO) systems in reference 5. 

A region of ultimate boundedness is defined as a compact region containing the 
origin in the state space such that trajectories starting from any initial state 
enter that region within a finite time and remain inside that region thereafter 
(ref. 6). When the nonlinearity in an SISO Lur'e-type system satisfies the sector 
condition in regions away from the origin, but violates it in a bounded region con- 
taining the origin, it was proved in reference 7 that there exists a region of ulti- 
mate boundedness, and an estimate of the region was obtained. 

The problem considered in this paper is somewhat similar to the Lur'e problem 
for MIMO systems. It differs from the Lur'e problem because the optimal feedback 
gain is included in the loop and because optimal LQ regulators have certain special 
properties. In particular, since the Riccati matrix is positive definite (under 
certain stabilizability and detectability assumptions (ref. 1)), its quadratic form 
represents a natural candidate for the Lyapunov function. By capitalizing on the 
special properties of optimal LQ regulators, estimates are obtained for the regions 
of attraction and ultimate boundedness for systems having nonlinearities belonging 
to the classes that were described previously, referred to as cases (1) and (2) 
hereafter. 

The organization of the paper is as follows. The formulation of the problem is 
given in the next section. Then, for case (1 ), a theorem is proved which gives the 
region of attraction followed by a useful corollary that is applicable when the con- 
trol weighting matrix is diagonal. The subsequent section includes theorems which 
give the regions of ultimate boundedness for case (2). The analytical results are 
demonstrated by application to attitude control of a large space antenna. On the 
basis of the analytical and numerical results, it is concluded that LQ regulators 
with better tolerance to nonlinearities can be designed by adjusting the performance 
function parameters. 
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point of intersection of hysteresis nonlinearity with a-axis 
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scalar defined in equation (18) 
regulator gain matrix 
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scalar defined in equation (37) 
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order of control vector 
m-vector valued function 
system order 
Riccati matrix 
state weighting matrix 
matrix defined in equation (32) 
modal amplitude for ith mode 
control weighting matrix 
space of real m-tuples 

entry in ith column of diagonal matrix R 
inverse images of E a and 
set defined in equation (36) 

sets defined in equations (17), (28), and (38) 
actual control vector 
command control vector 
Lyapunov function 

scalar proportional to volume of a hyperellipsoid 
axis system for hoop-column antenna (fig. 4) 
state vector 

ith component of state vector x 

maximum angular displacement about Y-axis of hoop-column antenna 

degree of stability 

matrix satisfying equation (6) 

positive scalar 

rigid-body pitch angle 

ith eigenvalue of P 
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Notation: 


smallest and largest eigenvalues 
scalar defined in equation (35) 

Lagrange multiplier 

damping ratio for the ith mode 

sets in Euclidean space of real m-tuples 

m-vector argument of N 

ith component of vector a 

nonlinear function defined in equation (11) 
ith component of vector <J> 

mode slope for the jth mode at location of ith actuator 
sets defined in equations (32) and (33) 
natural frequency of ith structural mode 


U 

n 

C 

€ 

[1 ,m] 


det( ) 
II II 


union of sets 
intersection of sets 
is a subset of 
is an element of 

the set of integers from 1 to m including 1 and m; e.g., 
means 1 < i < m 

determinant of a matrix 

Euclidean norm 


i C [1 ,m] 


Superscripts T, -1, and c respectively denote matrix transpose, inverse, and 
complement, A bar over a symbol denotes the boundary of a set. A dot over a symbol 
denotes the derivative with respect to time. 


PROBLEM FORMULATION 


The system is given by 


x = Ax + Bu 


(1 ) 


where x and u are n- and m-dimensional state and control vectors, and A and 
B are n x n and n x m constant matrices. Unlike the Lur’e problem (ref. 3), A 
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need not be Hurwitz; that is , A may have eigenvalues with nonnegative real parts 
Consider the infinite duration regulator problem where the following performance 
function is minimized: 


, r 2a t T T . 

J = J e (xQx+uRu) 

0 


dt 


( 2 ) 


where a is a nonnegative scalar representing the required degree of stability, 

Q is an n x n symmetric, positive semidefinite matrix, and R is an m x m positive 
matrix. The control vector u(t) which minimizes J in equation (2) is given by 
(ref. 1) 


u = Gx 


(3) 


where 


-1 T 

G = -R BP 


(4) 


and 


A T P + PA + 2aP + Q - pbr -1 b T p = 0 


(5) 


Since Q is positive semidefinite, it can be expressed as 


T 

Q = r r 


( 6 ) 


where T is an n x n matrix. Riccati equation (5) has a unique positive definite 
solution P if (A, B) is controllable and (T,A) is observable. These controllability 
and observability conditions are assumed to be satisified for systems considered in 
this paper. Under these conditions, the eigenvalues of (A + BG) have real parts less 
than -a • Let 


T 

V(x) = x Px 


(7) 


It can be shown that (ref. 1) 


v[ x( t ) ] < e 2(Xt v[ x( 0) ] 


( 8 ) 


That is, the closed-loop system has the degree of stability a. 
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( 9 ) 


In practical situations, nonlinearities exist in control actuators. In that 
case, equation (3) is replaced by 


u = Gx 
c 


u = N(u ) 
c 


( 10 ) 


where u c and u represent the commanded and actual control inputs, and N(a) 
denotes an m-vector valued, possibly time-varying nonlinear gain function of the 
m-vector argument a. From references 1 and 2, the closed-loop system is asymp- 
totically stable in the large (ASIL) if a T R[N(a) - 0.5a] > 0; that is , if 
a T R (f)(a) > 0, where 


<Ka) = N(a) - 0.5a (n ) 


The closed-loop system is given by 


x = x + B 4> (u c ) 


( 12 ) 


where u c is given by equation (9) and 


A = A + ( 1 / 2 ) BG 


(13) 


It has been established in reference 1 that A 1 is a strictly Hurwitz matrix. 


REGIONS OF ATTRACTION 

Consider nonlinearities that belong to the (0.5,®) sector in at least a neigh- 
borhood of the origin (case 1 described previously); that is, suppose the nonlinear- 
ity <]> (a) is such that the condition 


T 

a R 4> ( a ) > 0 (14) 


is satisfied for aCE a CR m , where £ a denotes a nonempty region containing the 
origin in the space R m of real m- tuples. For this case, an estimate of the region 
of attraction is obtained in this section. 

Suppose the region £ a contains a neighborhood of the origin (that is, the set 
{z|z G R m , II zll < 6} for some 6 > 0, where II zfl denotes the Euclidean norm). The 
inverse image S-j C R n of E a is defined as 
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S = {x | Gx G E } 
1 1 1 a J 


( 15 ) 



Let E denote the boundary of E-, , and let s be its inverse image; that is s„ 
a a 1 '1 

is the boundary of S -j : 

S = {x|Gx C E } (16) 

1 1 1 a J 


The following theorem gives an estimate of the region of attraction. 

Theorem 1 : If condition (14) is satisfied for 0 G E a , the closed-loop system 

of equations (1), (4), (9), and (10) is asymptotically stable (AS), and S a is a 

region of attraction where 


s a = {x|x t px < a} 


(17) 


T 

d = min (x Px) 


xGS 


1 


(18) 


Furthermore, the system has the degree of stability a inside S a . 

Proof : Differentiating V(x) in equation (7) and using equations (12), (13), 
(4), (5), and (9) results in 


V = -x T (Q + 2aP)x - 2 u T R <J>(u ) (19) 

c c 

Condition (14) is satisfied for u c £ E a , that is, for x G S-j where S-j is given 
by equation (15). Since E a contains a neighborhood of the origin of R m , S-| 
contains a neighborhood of the origin of R n . The situation under consideration is 

depicted in figure 2 for a two-dimensional system. Consider the region 

R§ = {x|v(x) < 6} . If R§ C S-| , R§ is a region of attraction because V > 0 
and V < 0 along all trajectories in R§. Since S^ contains a neighborhood of 
the origin, there exists a 6 > 0 such that R§ C S-\ • If 6 is successively 
increased, the largest 6 for which R§ C S 1 occurs when a boundary point of S-| 

is reached. This boundary point is also the value of x that minimizes V(x) for 

x G S v as stated in equation (18). Thus S a is a region of attraction. Also, from 
equation (19), V < -2aV inside S a , which implies that (ref. 1) 

v[ x( t ) ] < e v[x(0)] 


Thus the closed-loop system has the degree of stability a inside s a . 

In practice, the case when each component of N (and therefore, each ) 

is a function only of is more meaningful. Instead of condition (14), suppose 

the nonlinearities satisfy the condition 

a &> (a ) > 0 For I . < a. < i C [1,m] (20) 

111 1 1 1 2l 
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For example, the input-output graph of the saturating amplifier of figure 1(a) satis- 
fies condition (20). When condition (20) is satisfied, £ a is a region bounded by 
hyperplanes in R m given by 

S a = ^°^1i < a i < A 2i' 1 C (21) 

where i 1;L < 0 and l 2 i > °* Condition (14) is satisfied for a £ £ a . Let £ 
denote the set a 


£ a = {o|cr £ £ a , CL = or X^, i C [1,m]} (22) 

Corollary 1*1 : Suppose that R is a diagonal matrix with entries r^ > 0 
(i £ [1,m]) and the nonlinearities are such that condition (20) is satisfied. Then 
an estimate of the region of attraction for the closed-loop system given by equa- 
tions (1), (4), (9), and (10) is given by S a where 


S = {x|x T Px < d} 
a 

d= min [ ( X r .) Vb^Pb ] (23) 

iGt 1 ,m] 31 1 ' 1 1 

j£(1,2] 


where b ^ denotes the ith column of B. The system has the degree of stability a 
inside S a . 

Proof : In this case, differentiating V(x) in equation (7) results in 


V = -x T (Q + 2aP)x 


m 

l 

i=1 


r. u . <J>. (u . ) 
l ci 1 ci 


(24) 


where 


u . 
ci 


T 

g i x 


(25) 


and g^ is the ith row of G, given by 


T -IT 

< 3 i = ~ r i b ± P (26) 
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T 

The set s in this case consists of portions of hyperplanes g.x = JL 


nx 


(j C [1,2], i £ [1,m]), and is the region which is partially bounded by 

these hyperplanes# Using theorem A1 given in the appendix, 


min 

T 

g . x=Jl . . 


T 

(xPx) 



g i 


(27) 


Substituting equation (26) and performing minimization over the region S-| yields 
equation (23) in the statement of the corollary. As was seen in the proof of 
theorem 1, V > 0 and V < 0 along all trajectories in S a (because of assumed 
observability of (T, A)), which is defined by equations (17) and (23); thus, S a is 
a region of attraction# It is straightforward to see from equation (24) that since 
u ci ^i ( u ci ) > 0 in S a , the system has the degree of stability a inside S a # 

Corollary 1 #1 enables one to readily determine an estimate of the region of 
attraction for an LQ design, given Furthermore, equation (23) provides a 

method of adjusting weights in order to make regions of attraction larger. For 
example, if is small compared with the other (i.e., if the kth nonlinear- 

ity violates the sector condition much closer to the origin than the other nonlinear- 
ities do), one may increase the weight r^ to make the region of attraction larger. 


REGIONS OF ULTIMATE BOUNDEDNESS 

This section considers nonlinearities N(a) that lie outside the (0.5,®) sector 
only in a neighborhood of the origin (case 2 described previously). 

Let Eb C R m denote a compact region containing the origin. Suppose the func- 
tion N(a) is such that condition (14) is satisfied in z£, where the superscript c 
denotes the complement. Suppose that $(o) is bounded for a £ Let S 2 C R m 

denote the inverse image of that is, 

s 2 = {x|gx g I b } 


and let S 2 denote the boundary of It is assumed in this section that Q is 

chosen to be positive definite if a =0. The following theorem gives an estimate of 
the region of ulitmate boundedness. 

Theorem 2 : If condition (14) is satisfied for a and if 4> ( cr ) is bounded 

in Zb* then the region Sb is a region of ultimate boundedness for the closed-loop 
system given by equations (1), (4), (9), and (10), where 

S b = {x|x T Px < h} (28) 
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where 


T, 


h = max (x Px) 

(29) 

Q 1 UQ 2 


Q-) = {x|x T Qx = h, x C S 2 } 

(30) 

Q 2 = {x|x t Qx < )i, x C S 2 } 

(31) 

Q = Q + 2aP 

(32) 

H = -2 min [a T R <J)(a)] 

(33) 


oa k 

b 


Proofj As in the proof of theorem 1, differentiating V(x) in equation (7) 
results in 


V = -x t Qx - 2u£r <f> (u c ) ( 34 ) 

where Q is defined in equation (32). Since (by assumption) Q is chosen to be 
positive definite if a = 0, Q is positive definite. Then, condition (14) is 
satisfied for o c G££, that is, for x CSj. Therefore, from equation (34), 

V < 0 for x C S 2 . Using equation (33) in equation (34), we have 


V < -x T Qx + p, ( 35 ) 


Therefore, V < 0 for x C S 3 , where 


S 3 = {x|x t Qx > [i} 


(36) 


Therefore, V < 0 along all trajectories for x G S 3 U According to refer- 

ence 6 , the system is ultimately bounded in a compact region containing the origin 
of the form 

S b = (x| V(x) < h} 

• _ 

if V > 0 and V < 0 along all trajectories in s]~, and if v(x) -* » as tlx II -> a.. 
The least conservative estimate of the region of ultimate boundedness is obtained by 
finding the smallest h that satisfies these conditions. Thus the smallest hyper- 
ellipsoid containing the region (S 3 U S 2 ) c = S 2 D s| must be found. The smallest 
hyperellipsoid containing a region is the one containing its boundary, which in this 
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case is Q-j U Q 2 where and Q 2 are defined in equations (30) and (31). The 

situation is shown in figure 3 for a two-dimensional system. 

The maximization of equation (29) is difficult to perform analytically because 
of the nature of the region Q-j U £22* A more easily obtainable, but possibly more 
conservative (i.e., larger), estimate of the region of ultimate boundedness is 
obtained by performing the maximization over a larger region, that is, by using h* 
instead of h where 


h' = max (x T Px) 

rp A 

x i Qx=ji 


(37) 


Theorem 3 ; For the closed-loop system of equations (1), (4), (9), and (10), an 
estimate of the region of ultimate boundedness is given by 

= {x|x T Px 4 h ' } (38) 

where 

h' = \x/[2a + \ m (P -1 Q)] (39) 

where \x is defined in equation (33) and ^ denotes the smallest eigenvalue. 

Proof : The proof is obtained by minimizing V(x) over the region 

{x|x T Qx = |i}, which contains the region Q ^ U Q 2 used in theorem 2. Using 
theorem A2 in the appendix, 

max (x T Px) = n ^(Q" 1 ?) = p/X m (p- 1 Q) (40) 

x T Qx=)i 


where and X m denote the maximum and the minimum eigenvalues. (Since Q > 0 

and P > 0, the eigenvalues of Q -1 P are all real and positive (ref. 8).) Equa- 
tion (39) is obtained by using equation (32) in equation (40). 

When R is diagonal, we have from equation (33) 


m 

\i = -2 min £ 

aCZ, i=1 
b 


r. a. 4>. (a. ) 

iii i 


(41) 


For making or small, \i (a positive scalar) should be made as small as 

possible. From equation (41), \i can be made smaller by reducing the weights r^ 
corresponding to those input channels having nonlinearities that most severely 
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violate the (0.5,°°) sector condition (i.e., for most negative in 

the violation region). Also, from equation (39), h' can be made smaller by 
(1) increasing a or (2) increasing Q. Although the reduction of is not 

guaranteed because of the dependence of P on Q and <x, these procedures provide 
methods of performing parametric studies to obtain an LQ design with an acceptably 
small region of ultimate boundedness. 


NUMERICAL RESULTS 

In order to demonstrate the analytical results obtained, the problem of control- 
ling the rigid-body pitch angle and elastic motion of a large, flexible space antenna 
was considered. Figure 4 shows the hoop-column antenna concept, which basically con- 
sists of a deployable central mast attached to a deployable hoop by cables held in 
tension. A secondary drawing surface is formed by quartz or graphite stringers 
attached between the hoop and the mast, and the radio-frequency (RF) reflective mesh 
is attached to the secondary drawing surface by mesh shaping ties. To achieve the 
required RF performance, the rigid-body attitude of the antenna must be precisely 
controlled, and the elastic motion must be kept very small. The mathematical model 
considered in this paper includes the rigid-body pitch angle about the Y-axis and the 
first two bending modes in the XZ-plane. The nominal pitch-axis model is given by 


10 


T, + To 


(42) 


where I is the Y-axis moment of inertia, 0 is the rigid-body pitch angle, and 
T 1 and T 2 are the Y-axis control torques applied by control moment gyros (CMG's) 

points 1 and 2 shown in figure 5. The elastic motion for the ith structural mode 
is given by 


qi. + 2 Pi W i<Ji + W iqi = ^ r 1i T 1 + ¥ 2 .T 2 


(43) 


where q if p^, and denote the modal amplitude, inherent damping ratio, and the 

natural frequency for the ith mode, and Yji denotes the ith mode slope at actuator 
location j. The parameters of the 1 22-m-diameter hoop-column antenna are taken from 
reference 9 and are given in table I. The elastic deformations due to the two bend- 
ing modes are shown in figure 6. An optimal LQ regulator can be designed for this 
model to minimize the performance function in equation (2). It is assumed in this 
example that the entire state vector x, defined as 


• • • T 

X = (0,0,q i ,q i ,q 2 ,q 2 ) (44) 

is available for feedback. The following two types of nonlinearities are considered: 
(1) saturating actuators and (2) actuators with hysteresis. 
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Saturating Actuators 


Assume that the actuator characteristics of both actuators are as shown in fig- 
ure 1(a) and are given by 


T 


c 


sgn(T 


)T 

max 


For 

For 



< T 

max 

> T 
max 


(45) 


where 


c r 


T, and T, 


max 


denote the command torque, the actual torque, and the 


maximum torque, and sgn( ) denotes the signum function. The value of T max 
is assumed to be 1.627 x 10^ N-m. (This value is rather large from a practical 
viewpoint, but is only used as an example.) Since the slope of the nonlinearity is 


> 2T 


max* 


unity in the linear region, the (0.5,®) sector is violated for |l 
nominal LQ regulator was first designed to obtain a rigid-body closed-loop frequency 
and damping ratio of 0.138 rad/sec and 0.707, respectively, and closed-loop struc- 
tural mode damping ratios of at least 0.5. The weighting matrices and the closed- 
loop eigenvalues for the nominal design are shown in table II. The degree of 
stability parameter a was assumed to be zero for the nominal design. In the pres- 
ence of control saturation, there exists a bounded region of attraction S a . The 
larger the S a , the better the tolerance of the design to the nonlinearities. Since 
the estimate of S a is a hyperellipsoid, it is difficult to visualize S a . One 
measure of the size of S a would be the maximum value of the angular displacement 
(rigid-body plus elastic) at one of the sensor locations. Suppose an attitude sensor 
is placed at the location of actuator 1 • The angular displacement about the Y-axis 
is then given by 


T 

y = c x 


(46) 


where the constant vector is 


T 

c 


(bO^^O^^O) 


(47) 


Therefore, according to theorem A3 of the appendix, the maximum angular displacement 
y m ax within S a is 

y = \1c T P cd (48) 

max 

Another measure of the size of S a is its volume. The volume of a hyperellip- 
soid is proportional to the product of its semimajor axes. For the hyperellipsoid 
given by 


T 

x Px = d 


( 49 ) 
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it is sh own in theorem A4 of the appendix that the ith semimajor axis is given by 
\|d/\^ ( P ) , where \^(P) denotes the ith eigenvalue of P. Thus the volume is pro- 
portional to V' where 


V 1 



n x. (p) 

i 



(50) 


where det ( ) denotes the determinant of a matrix. These two measures, y m ax and 
V', are used herein to evaluate the LQ designs. 

To investigate the effect of variation of R on the size of S a , a series of 
LQ regulators were designed, with R increased by a factor of \JT6 at each step. The 
initial value of R was 0.01 times its nominal value (identity matrix). Figure 7 
shows a plot of y max and V 1 as R increases. For the nominal design, y max 
is 68.5°, which is satisfactorily large. As R increases, both y max and V 1 
increase, because increasing R decreases the magnitude of control effort, which is 
thus less likely to reach the saturation limits. Of course the performance generally 
deteriorates as R increases. 


Now assume that the two saturating actuators have different saturation limits. 
The saturation limit for actuator 1 is T max as before, but that for actuator 2 is 
0.125T max . For this case, the nominal LQ design yielded y max of only 8.6°. As 
discussed previously, equation (23) suggests that increasing the control weight r 2 
for actuator 2 should increase the size of S a . Therefore, r 2 was increased by a 
factor of 2 at each step and a series of LQ regulators were designed. As shown in 
figure 8, both Y max and V' increase as r 2 increases up to the sixth step. 

To investigate the effect of the degree of stability parameter a on S a , a 
was increased by 0.02 at each step for 20 steps, and y max and V 1 were computed 
for the resulting LQ designs. As shown in figure 9, both y max and V* decrease 
with increasing a. Because increasing a increases feedback gains, the control 
effort reaches the saturation limits earlier. 


Actuators With Hysteresis 

Now assume that the actuator characteristics for each torque actuator are as 
shown in figure 1(c), with K h = 1 and a = 0.25 N-m. For this case, \i was deter- 
mined from equation (41) as | x - a 2 (r-j + r 2 ) • The region of ultimate boundedness 

is given by equations (38) and (39). The smaller the size of S£, the better the 
tolerance of the design to these nonlinearities. 

The nominal LQ design described in the previous section yielded y max = 2.32° 
for the region of ultimate boundedness. As in the previous section a series of LQ 
regulators were designed by increasing R by a factor of \JTo at each step, and 
y m ax and v ' were computed for the region of ultimate boundedness for each LQ 
design. As shown in figure 10, both y max and V' increase as R increases. 

To investigate the effect of a on S^, a was increased at each step by 0.02 
(starting with a = 0), and y max and V* were computed for each resulting LQ 
design. Figure 11 shows the variation of y max and V' with a. There is a large 
decrease in both y max and V* from the first to the second step, because 
[ 2a + \ m (P -1 Q)] appears in the demoninator in equation (39) and a is large com- 
pared with X m (P“ 1 Q). As a increases further, y max and V' both continue to 
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decrease, but at a slower rate. From this example, a appears to be an important 
design parameter, which can significantly affect the size of the estimate of the 
region of ultimate boundedness. 


CONCLUDING REMARKS 

Closed-loop stability was investigated for multivariable linear time-invariant 
systems controlled by optimal LQ regulators when nonlinear gains are present in the 
control channels. Two types of nonlinearities were considered: (1) nonlinearities 

that lie in the (0.5,®) sector at least in a neighborhood of the origin and (2) non- 
linearities that lie in the (0.5,®) sector in regions away from the origin, but per- 
haps escape that sector in a neighborhood of the origin. Making use of the special 
properties of LQ regulators, estimates were obtained for the region of attraction for 
case (1) and for the region of ultimate boundedness for case (2). The sizes of these 
regions represent measures of robustness of a given LQ regulator when particular 
types of nonlinearities are present in the feedback loop. The expressions obtained 
also provide methods for selecting the performance function parameters (i.e., the 
state and control weighting matrices and the degree of stability) in order to design 
LQ regulators with better tolerance to nonlinearities, that is, to obtain a larger 
region of attraction or a smaller region of ultimate boundedness. The analytical 
results obtained were demonstrated by application to the problem of controlling the 
pitch-axis attitude and elastic motion of a large, flexible space antenna. Based on 
the analytical and numerical results, it was concluded that 

(1) Decreasing the state weighting matrix Q (or, equivalently, increasing the 

control weighting matrix R) results in larger region of attraction and 
larger region of ultimate boundedness. 

(2) Increasing the degree of stability parameter a decreases the size of both 

the region of attraction and the region of ultimate boundedness. In 
particular, by choosing a to be a small positive scalar instead of zero, 
the region of ultimate boundedness can be made significantly smaller. 

(3) For case (1), if a nonlinearity escapes the (0.5,®) sector closer to the 

origin than the other nonlinearities, the region of attraction can be made 
larger by increasing the weight corresponding to that control channel. 

This paper assumed that the complete state vector is available for feedback. In 
practice, however, only a few sensor outputs (fewer than the dimension of the state 
vector) are usually available, and an observer or a state estimator must be used. 
Future research efforts should be directed toward the synthesis of robust overall 
controllers, which include an LQ regulator and a state estimator. It is expected 
that the analytical results obtained in this paper would be useful toward achieving 
that goal. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 3, 1984 
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APPENDIX 


DERIVATION OF ALGEBRAIC OPTIMIZATION RESULTS 

This appendix contains the derivation of some algebraic optimization results 
used in this paper, as well as the expression for the semimajor axes of a hyperellip- 
soid. Throughout the appendix, x and c denote real n-vectors, P and Q denote 
real symmetric n x n positive definite matrices, and A, \i, and d denote real 
positive scalars. 

Theorem A1 : 


T 2 T -1 

min (x Px) = A /c P c (A1 ) 

T o 

c x=A 


Proof : With the Lagrange multiplier denoted by v, the Hamiltonion is given by 

T T 

H = x Px + v(A - c x) 

Therefore, the necessary condition for minimum is 

5H/5x = 2Px - cv = 0 (A2) 

Premultiplying the above equation by x T , using the constraint (c T x = A) , and solving 
equation (A2) for v results in 

v = 2J 1 /A (A3 ) 

where J 1 = x T Px. Substituting equation (A3) in equation (A2) and solving for x 
yields 

x = P c^/A ( A4 ) 

Therefore , 

T 2 T -1 2 

J, = x Px = P c/A (AS) 

from which equation (A1 ) is obtained. 
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Theorem A2 : 

T -1 

max (x Px) = \ (Q P) (Jt 
M 

T^_ 

x Qx=|i 

where ) denotes the largest eigenvalue. 

Proof: The Hamiltonion is given by 

T T 

H = x Px + v([l - x Qx) 

Therefore, the necessary condition is 

9H/9x = 2Px - 2vQx = 0 

Equation (A8) can be rewritten as 

Q 1 Px = vx 

Therefore, v is an eigenvalue of Q -1 P. 
using the constraint (x T Qx = ji) yields 

J 1 = X ( Q _1 P ) \i 

Therefore the maximum value of J-j is given 
Theorem A3 : 

, T . \|~T — 1 " 

max (c x) = Nc P cd 

T 

x Px=d 

Proof: The Hamiltonion is given by 

T T 

H = cx+v(d-xPx) 

The necessary condition is 

9H/dx = c - 2vPx = 0 


(A6) 


(A7) 


(A8) 


(A9) 

Premultiplying equation (A8) by x T , and 


(A10) 

equation (A6). 


(All) 


(A12) 


(A13) 
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Premultiplying equation (A13) by x T and using the constraint (x T Px = d) yields 


v = J 2 /2d (A14) 

where J 2 = c T x. Substituting equation (A14) in equation (A13) and solving for x 
results in 

x = P _1 c(d/J 2 ) (A1 5 ) 

From equation (A15), 

T T -1 

J 2 = c x = c P c(d/J 2 ) (A16) 

Equation (All) is obtained from equation (A16). 

Theorem A4: For the hyperellipsoid given by 

T 

x Px = d (A17) 


the semimajor axes are given by Njd/kjJP) , where ^(P) is the ith eigenvalue of 
P (i C [1,n]). 

Proof : Since P is symmetric, it is orthogonally similar to a diagonal matrix 
(ref, 8). That is, there exists a real orthogonal n x n matrix E such that 

T 

E PE = A p (A18) 

where Ap is the diagonal matrix with the eigenvalues of P as its entries. Let 
x = Ey, where y is an n-vector. Then equation (A17) becomes 

T. 

y A p y = d (A19) 

The ith semimajor axis is obtained by making yj = 0 for j * i and is given by 


= d/\.(P) 


(A20 ) 


which proves the theorem. 
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TABLE I.- HOOP-COLUMN ANTENNA PARAMETERS 


Inertia about Y-axis , 


kg-m 2 


First bending mode parameters: 

Pi 

, rad/sec 


Second bending mode parameters : 

P2 

0 ) 2 / rad/sec 


5.748 x 10 6 


0.01 

1.35 

-3.385 x 10 -3 
2.938 x 10 -4 


0.01 

5.78 

-1 .754 x ID' 5 
2.270 x 10~ 4 


TABLE II.- WEIGHTING MATRICES AND EIGENVALUES (NOMINAL DESIGN) 
FOR LQ REGULATOR FOR HOOP-COLUMN ANTENNA 


Weighting matrices: 


Q = diag (1.2 x 10 9 , 1.2 x 10 9 , 1.0 x 10 3 , 3.0 x 10 5 , 1.0 x 10 3 , 1.0 x 1 0 9 ) 
R = diag( 1 ,1 ) 


Closed-loop eigenvalues ( j = \PT ) : 


-0.0761 ± j 0.0759 
-0.800 ± jl .09 

-3.09 ± j 4 .88 


Feedback gain matrix: 


X 

CM 

. 

CM 

1 

10 5 

1 

u> 

• 

U) 

X 10 5 

3.31 

X 

o 

5.36 

x 10 2 

3.41 

X 

o 

(jO 

2.05 x 10 3 

-2.64 x 

10 5 

-3.53 

in 

o 

X 

5.89 

*"o 

X 

-2.65 

x 10 1 

r- 

CO 

. 

CM 

X 10 3 

-3.1 x 10 4 
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Figure 4.- Hoop-column antenna concept. 



Figure 5,- Assumed actuator locations on hoop-column antenna. 
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Figure 6*- Bending mode shapes of hoop-column antenna. 
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Figure 11.- Effect of increasing a on region of ultimate boundedness 
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